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Abstract 

This paper reports the tenth-order contributions to the g — 2 of the electron a e and those 
of the muon from the gauge-invariant Set 11(c), which consists of 36 Feynman diagrams, 
and Set 11(d), which consists of 180 Feynman diagrams. Both sets are obtained by insertion 
of sixth-order vacuum-polarization diagrams in the fourth-order anomalous magnetic moment. 
The mass-independent contributions from Set 11(c) and Set 11(d) are —0.116 489 (32)(a/7r) 5 and 
—0.243 00 (29)(a/7r) 5 , respectively. The leading contributions to an, which involve electron loops 
only, are -3.888 27 (90)(a/vr) 5 and 0.497 2 (65)(a/vr) 5 for Set 11(c) and Set 11(d), respectively. The 
total contributions of the electron, muon, and tau-lepton loops to a e are -0.116 874 (32)(a/vr) 5 for 
the Set 11(c), and —0.243 10 (29)(a/7r) 5 for the Set 11(d), respectively. The contributions of electron, 
muon, and tau-lepton loop to are —5.559 4(ll)(a/7r) 5 for the Set 11(c) and 0.246 5(65)(a/7r) 5 
for the Set 11(d), respectively. 

PACS numbers: 13.40.Em,14.60.Cd,12.20.Ds,06.20.Jr 
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I. INTRODUCTION 



The anomalous magnetic moment g—2 of the electron has played the central role in testing 
the validity of quantum electrodynamics (QED) as well as the standard model. The latest 
measurement of a e = (g—2)/2 by the Harvard group has reached the precision of 0.24 x 10 -9 



[i 



l2|: 



a e (HV08) = 1 159 652 180.73 (0.28) x 10~ 12 [0.24ppb] . (1) 
: theory consists of QED corrections of up to the eighth order 



At present the best prediction o 

nil nn 

[3H5J, and hadronic corrections [6Hl2| and electro- weak corrections [13Hl5l| scaled down from 
their contributions to the muon g — 2. To compare the theoretical prediction with the 
experiment (JTJ) , we also need the value of the fine structure constant a determined by a 
method independent of g—2 . The best value of such an a has been obtained recently from 
the measurement of h/m-^, the ratio of the Planck constant and the mass of Rb atom, 



combined with the very precisely known Rydberg constant and m^K/m, |16l|: 

or -1 (RblO) = 137.035 999 037 (91) [0.66ppb]. (2) 

With this a the theoretical prediction of a e becomes 

a e (theory) = 1 159 652 181.13 (0.11)(0.37)(0.77) x 10~ 12 , (3) 

where the first, second, and third uncertainties come from the calculated eighth-order QED 
term, the tenth-order estimate, and the fine structure constant (j2J), respectively. The theory 
([3]) is thus in good agreement with the experiment ([1]): 

a e (HV08) - a e (theory) = -0.40 (0.88) x 10" 12 , (4) 

proving that QED (standard model) is in good shape even at this very high precision. 

An alternative test of QED is to compare the a of fl2]) with the value of a determined 
from the experiment and theory of g — 2 : 

a -\a e 08) = 137.035 999 085 (12)(37)(33) [0.37ppb], (5) 

where the first, second, and third uncertainties come from the eighth-order QED term, 
the tenth-order estimate, and the measurement of a e (HV08), respectively. Although the 



uncertainty of a _1 (a e 08) in (J5J) is a factor 2 smaller than a _1 (RblO), it is not a firm factor 
since it depends on the estimate of the tenth-order term, which is only a crude guess [17] . In 
anticipating of this challenge we launched a systematic program several years ago to evaluate 



the complete tenth-order term 



18- 



The tenth-order QED contribution to the anomalous magnetic moment of an electron 
can be written as 



,(io) 



«N 5 
7T 



+ A { 2°\m e /m^) + A { ™ ] \m e /m T ) + A^ 0) (m e /m tlJ m e /m 7 



(6) 



'he result of 105 



25|]. Five of the 



where the electron-muon mass ratio m e /m M is 4.836 331 71 (12) x 10 -3 and the electron-tau 
mass ratio m e /m T is 2.875 64 (47) x 10 -4 17J. The contribution to the mass-independent 
term A^ coming from 12672 Feynman diagrams may be classified into six gauge-invariant 
sets, further divided into 32 gauge-invariant subsets depending on the nature of closed lepton 
loop subdiagrams. Thus far, results of numerical evaluation of 24 gaug e-invariant subsets, 
which consist of 2785 vertex diagrams, have been published 0, 1211 24]. 
vertex diagrams of Set I(i) has been recently submitted for publication 
subsets had also been calculated analytically 2^1 • Our calculation is in good agreement 
with these analytic results. 

In this article we report the evaluation of contributions of two gauge-invariant subsets, 
Set 11(c) and Set 11(d), which consist of fourth-order vertex diagrams containing vacuum- 
polarization subdiagrams of sixth order. The effect of insertion of a gauge-invariant set of 
closed lepton loops in an internal photon line of momentum q is expressed by the renormal- 
ized vacuum-polarization tensor of the form 



n^(?) = (9Y-g 2 r)n(g 2 ) 



(7) 



where the scalar vacuum-polarization function II(g 2 ) vanishes at q 2 = on carrying out the 
charge renormalization. 

The Set 11(c) consists of 36 Feynman diagrams. A typical diagram of this set is shown 
on the left-hand side of Fig. [TJ It is obtained by insertion of proper sixth-order vacuum- 
polarization diagrams containing two closed lepton loops (see Fig. [2]) in the fourth-order 
anomalous magnetic moment M^ a or My, represented by Fig. [3j These diagrams can be 
represented by 4 independent integrals taking account of various symmetry properties. 

The Set 11(d) consists of 180 Feynman diagrams. A typical diagram of this set is shown 




11(c) 



11(d) 



FIG. 1: Typical diagrams of the tenth-order Set 11(c) and Set 11(d). 





4a(P2) 



4b(P2) 



FIG. 2: Sixth-order vacuum-polarization diagrams consisting of two fermion loops. 



on the right-hand side of Fig. [TJ It is obtained by insertion of proper sixth-order vacuum- 
polarization diagrams consisting of one closed lepton loop in the fourth-order anomalous 
magnetic moment. These diagrams can be represented by 16 independent integrals taking 
account of various symmetry properties. 

Evaluation of the contribution of the Set 11(c) to the mass-independent term is 
straightforward since an exact spectral function II^ 4 ' 2 ) for the diagrams of Fig. [2] is known 
for the diagrams whose two lepton loops have the same mass 28j]. However, evaluation of 
the mass- dependent term A 2 requires function of m e / or m e j m T , which is not 

available at present. In order to cover both cases, we follow an alternative approach 29( in 
which we construct a Feynman-parametric integral of the sixth-order vacuum-polarization 




FIG. 3: Self-energy-like diagrams of fourth order. 
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FIG. 4: Sixth-order vacuum-polarization diagrams consisting of a single fermion loop. 

function II^ 4 ' 2 ) and insert it in the virtual photon lines of the Feynman-parametric integral 
of the forth-order anomalous magnetic moment M4. 

For the Set 11(d) an exact vacuum-polarization function II^ (see Fig. H]) is not known, 



although the Pade approximant is known to provide a good approximation |30H32| . We 
follow here primarily the approach [29| which utilizes the vacuum-polarization function H^ 6 ) 
itself, instead of its spectral function. The calculation utilizing the Pade approximant of the 
spectral function is also carried out to provide an independent check. 

Parametric representations of several vacuum-polarization functions are presented in 
Sec. Ull where explicit definitions of functions are given. As an illustration of our approach, 
insertion of vacuum-polarization function in M 2 is presented in Sec. IHIl Insertion of U^ 2 \ 
n( 4 ) ; n^ 4 ' 2 ) and Il( 6 ) in M 4 is described in Sec. IIVI Although most results of Sec. [TT1 and 
Sec. IHIl are concerned with quantities of low orders, they are needed in carrying out the 
renormalization of the tenth-order terms. In cases where the spectral function is avail- 
able, we present alternative ways which provide a consistency check of the numerical work. 
Application of these methods to Set 11(c) and Set 11(d) is described in Sec. [V] and Sec. IVTl 
Summary and concluding remarks are presented in Sec. lVIII For simplicity the factor (a/V) 5 
is omitted in Sees. [TTI — PVTI 
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FIG. 5: Fourth-order vacuum-polarization diagrams consisting of one lepton loop. 

II. PARAMETRIC REPRESENTATION OF VACUUM-POLARIZATION FUNC- 
TION 



As is shown in Ref. 29] the second-order vacuum-polarization function can be written in 
the form 

n< 2 >M = /(<fe)^in(^), (8) 

with 



(dz) = dz 1 dz 2 S(l - zi 2 ), U = z 12 , V = zi 2 m\, V — V - xG, x = q 2 , 

G = z 1 A 1 , A 1 = z 2 /U, D = 2A 1 (1-A 1 ), (9) 

where z\ and z 2 are Feynman parameters of leptons forming the loop, Z\ 2 = zi + z 2 , and ni\ 
is the rest mass of the lepton. 

As a preparation for constructing II^ 4 ' 2 ) let us first construct the parametric integral of 
the fourth-order vacuum-polarization function il 1 - 4 ^. It has contributions from one diagram 
Pi a and two diagrams of Fig. [5j 

The contribution of P ia , renormalized at q = but with subvertex divergences not yet 
removed, can be written as |29j 



n (4a) . 



x 



{dz) 



£o(i_ n a (Vq 

U 2 \V V J f/ 3 n \V 



(10) 



where z%, z 2 , z$, Z4 are Feynman parameters for the electron lines and z a is that of the photon 
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line and 



(dz) = dzidz 2 dz 3 dz 4 dz a 5(l - 21234a), z% > 0, 

B\\ = Z23a, B12 = Z a , B 2 2 = Z 44a , U = Z\±B\\ + Z23-B12, 

Ai = (z 3 B 12 + z 4 B n )/U, A 2 = {z 3 B 22 + z 4 B 12 )/U, 

A 3 = A 2 - 1, A A = A l - 1, 

Vq = z 1234 ml + A 2 z a , G = Z\A\ + z 2 A 2) V = V - xG, x = q 2 , 

D = ((A 1 + A 4 )(A 2 +A 3 )-A 1 A 4 -A 2 A 3 )m 2 l , 

D x = (A 4 A 2 + A 3 A 4 )B l2 - A 1 A 4 B 22 - A 2 A 3 B U . 



(11) 



Here Z23a = z 2 + z z + z a, etc., and A is the (infinitesimal) photon mass. 

This integral contains ultraviolet (UV) divergences arising from the vertex diagrams { 
2, 3, a} and {1,4, a}, which can be removed by the i^-operation and ^14-operation 29] . 
respectively. The renormalized function iTftn. can be written as 



(12) 



where 



An (4a) = ^ _ R ^ _ ^ 14 ) n ( 4 «). 



(13) 



I[( 2 ) is the second-order vacuum-polarization function given by Eq. (jSJ), and = L 2 — L 



uv 
2 1 



where L^ v is the UV-divergent part of the second-order vertex renormalization constant L 2 
defined by the i^-operation and is the remainder including an infrared divergent part of 
L 2 . 



The renormalized vacuum-po 



similarly by the .fT-operation [29J. It leads to 



arization term coming from of Fig. [5] may be handled 



Uif^ = AU^-2BfU^\ 



where B^ = B 2 



5 2 uv [29| and 



An( 4b )( 



x) 



(dz)(l-K 2 ) 



Do(l_ 
U 2 \V 



1 \ xC D x fVo 

Vol u 2 v u 3 n \v 



(14) 



(15) 
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with 



B\l — Z2a, B\2 — Z a , B22 — ^134a, U — ^134-Bll + ^2-Bl2j 

A 1 = z 4 B u /U, A 2 = z 4 B 12 /U, A 3 = A U A i = A 1 - 1, 
G = z X3 Ai + z 2 A 2} V = V - xG, x = q 2 , 

D = (4Ai - A 2 )A 4 m 2 , C = -A 2 A 2 A 4 , D 1 = A i (A 1 + 3A 4 )B 12 . 



(16) 



Vo and (dz) have the same form as in pip . 

The function IT^ 4 ' 2 ) is obtained readily from by insertion of the spectral representation 
of the vacuum-polarization loop 



IT(g 2 



dt- 



p(t) 



(17) 



r Jo -<? 2 + 4m 2 /(l -t 2 y 

in the virtual photon line of il^ 4 -* carrying momentum q. This equation can be regarded as a 
superposition of vector particles of mass 4m 2 ,/ (1 — t 2 ), where m 2 is the mass of the inserted 
loop particle. The net effect is expressed as the replacement 



dt 



P(t) 



q 2 -4m 2 /(l -t 2 ) 



(18) 



For the second-order electron loop the spectral function is given by 29] 

t 2 (l-t 2 /3) 



P {2 \t) 



1-t 2 



(19) 



Do(]_ J_\ £i (Vo 

u 2 \v Vol u 3 n \v 



From (TlTJT) and (|T8|) we obtain 

An( 4a ' 2 )(x) = J dtp {2 \t) j (dz)(l - K 23 - K u 
where x = q 2 , and D and Di are given in Eq. ( TTTj) . Vo and V have the form 

V = V - xG, x = q 2 , V = z 12 3 4 ml + z, 
Similarly, from ( JT5|) and ffTH]) . we obtain 

An<«» W = /' / - ft) ( I - 1) + ^ + § i„ (| 

where D , C and are given in Eq. ffTB"]) . and V and Vq are given by ( 12~T|) . 



(20) 



4m 2 
1-t 



2 ' 



(21) 



(22) 
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In the same manner as in (fTUj) and (TP) 
n( 6 )(g 2 ) can be written in the general form 



the sixth-order vacuum-polarization function 
29] 



An (6) ( 



x 



(dz) 



u 2 V v 2 v 2 1 



JJ3 



1 

T 7 



i 



' 

xB\ 



x(B + xCq) 



Do 



U 2 V 2 



V'o 



f/ 3 V u 4 



~~\~ ~~t In I -7-7 



V 



where 



V'o — ^123456^1; 



(23) 



(24) 



and Z?o is m\ times Do, and Bo and Z)i are times Bo and Z?i, given in Ref. 33]. For 
simplicity i^-operation is not shown explicitly. For the sixth-order vacuum-polarization 
diagrams Pqc and Pqd of Fig. HI which contain 4th-order lepton self-energy subdiagrams, 
the i^-operation subtracts only UV-divergent parts Sm^ and Sm^ , which are different 
from the full mass-renormalization terms Sm^a and Sm^- This causes no problem for the 
renormalization of vacuum-polarization function which has no infrared divergence. However, 



the fully renormalized formula is simpler if the residual mass-renormalization term 5m 
Sni4 a — 5m^) and 6mf b (= Sm^ — Sm^ b y ) are also subtracted. This subtraction is performed 
by the /^-subtraction method introduced in Ref. 19j. We shall therefore include the R- 
subtraction operation in Eq. (123]) whenever it is needed. 



III. INSERTION OF VACUUM-POLARIZATION FUNCTION INTO THE 
SECOND-ORDER ANOMALY 

Before discussing insertion of vacuum-polarization (VP) diagram in M4, let us consider 
insertion in M2. The Feynman parametric representation of the second-order magnetic 



moment can be written in the form 



22| 



where y\ and y a are Feynman parameters of the electron line and the photon line, respec- 
tively, 

(dy) = dyxdyj^-yta), yi a = yi + y a , 
U y = y ia B n , B n = 1, A 1 = y a /U y , V y = yi- y x A x , F = y^l - A x ), (26) 
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and the rest mass of the open electron line (or, open lepton line) is chosen to be 1. The 
insertion of the vacuum-polarization function II into the photon line a can be written as 



M 



2,P 



dtp(t) 



(dy) 



m V y + Am\y a (l-V)^ 



(27) 



following Eq. (fTBj) . Comparing ([2~7]) with (fl7|) . we can write 



2,P 



dt- 



P(t) 



(dy) F 

Ul y a J ^V y /y a + Am\(l-t^ 
(dy) F II(g 2 ) 

U$ y a q 2 q 
(dy)F 



-Vy/ya 



U 2 y V y 



n(-V y /y a) 



(28) 



This gives us a simple and transparent recipe for insertion of the vacuum-polarization func- 
tion in the photon line a of M 2 : 



v y v y 

This is identical with Eq. (5.8) on page 285 of Ref. [29| noting that 

V 



(29) 



V = z 12 m x , V = V - xG\, 



x=-V y /y a i 



w 



(30) 



V-V 

Note that this derivation requires only the analytic property expressed by the spectral rep- 
resentation. No actual knowledge of the spectral function is required. 

Making use of Eq. ()29]) and the form of II (2) (x) given by Eq. (jSJ), M2 5 p 2 can be readily 
written in the form 



M- 



2,P 2 



(///(] - //) / (dz)^\n 



x = —- 



y a 



i-y 



(31) 



where y = y x , y a = 1 - y 1: and V = V - xG. 
In the case P = P4, we obtain from Eq. ([15 



M 



2,P 4 



dy(i -y) I (dz) 



D ( 1 1\ xC Di 

+ + — m 



(32) 



x=-y 2 /{l-y) 



_u 2 \v V J U 2 V u 3 

where V = Zi2u m \ an d V = V — xG. ^-operation (and i?-subtraction) are not shown 
explicitly for simplicity. This corresponds to Eq. (5.20) in page 288 of Ref. 29) noting that 

V 



V = z l2 m 11 V = V - xG\ x= _ Vy/ya , 
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W 



V-Vn 



(33) 



Similarly, for insertion of P = Pq, we have the general structure 



M, 



2,P 6 



dy(i -v) I (dz) 



U 2 \V 2 



+ 



Ei 

u 3 



i 

v 



1 



xBi D 2 , 



V 2 

Vo 
V 



+ 



xB + x 2 C 
U 2 V 2 



(34) 



x=-y 3 /(l-y) 

where Vq is given in Eq. (I24p . This equation, with = 1, corresponds to Eq. (5.32) on 



page 293 of Ref. 
diagrams. 



29] except that it includes .R-subtractions besides fT-operations for some 



IV. INSERTION OF VACUUM-POLARIZATION FUNCTION INTO A FOURTH- 
ORDER MAGNETIC MOMENT 



The fourth-order magnetic moment M4, which consists of two parts M^ a and M4&, has 
the form 



Ma 



(dy) 



E + Co N + Z + N 1 + Z 1 



U 2 V 2 
y v 



Wv 



(35) 



where (dy) = dyidy 2 dysdy a dybS(l — 2/12306), 2/1,2/2, 2/3 are Feynman parameters of lepton lines, 
2/ a , 2/b are Feynman parameters of photon lines, and 



V y = 2/123 + A 2 2/ab -G, G = y x A x + y 2 A 2 + y 3 A 3 , etc. 



(36) 



E ,Cq, N , Z , Ni, Zi are functions of Feynman parameters defined for M iq and M 46 , respec 



tively. Their explicit definitions are given in pages 266 and 267 of Ref. [29j. The tilde on 
Cq is introduced here to avoid confusion with Co introduced in the definition of II. When a 
VP function is inserted into the photon line a by using its spectral function representation, 
the denominator V y is replaced by V y + y a R(z). In the case of the second-order VP function 
P 2 , where z\,z 2 are Feynman parameters of two fermions forming the vacuum-polarization 
loop, we have R(z) = m\z\ 2 / (ziz 2 ). 

For terms of Eq. fl35|) proportional to 1/V y , we can apply the substitution rule fl29|) 
directly. For the term proportional to we may rewrite the denominator using the 
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formula 



(dz) 



d 



(V y + y a R(z)y 



dV, 



[dz] 



V y + y a R(z) 



d U(-V y /y a ) 

dVy Vy 



V 2 

v 



V„ 



This leads to the structure of the form 



4,P 



(dy) 



P + C iV + Z Ni + Z 



U^Vy 



?7 2 1/ 2 
3/ y 



+ 



U*Vy 



dV„ 



■n(-W) 



(37) 



iv + z an(-Vya) 



(38) 



For P = P2 we have 



M 4 , P2 = -]T /(dy) /(d 



Eq + Cq iV + Z 



£/ 2 \/ 2 
y y 




which follows from Eq. (jSJ) and 



1 

Va 



{dz 



DoG 
U 2 V : 



(39) 



(40) 



where Vq and V are given by Eq. ( 1331 . and 6 means the sum of insertions of P2 in photon 
lines a and b. Note that the structure of M 4 is largely kept intact by the insertion of II. 

For P = P 4 , where P 4 represents P 4a or P 4b , U^(—V y /y a ) is given by Eq. (}T5]) and its 
derivative can be written in the form 



dU^(-V y /y a ) 
dV„ 



Va 



(dz) 



DoG_ Co (}_ xG 

Ipv 2 + u 2 \v + v 2 



DiG 

Tpv 



(41) 



where, for P^, V — Vq — xG, Vq = z^^mf, G = zi 3 Ai + Z2A2, x = —V y /y a . Similarly for 
P = P 4o . Substituting Eqs. (fT5|) and (14T1) in Eq. (|38|) . we obtain M 4a (p 46 ), etc. To avoid 
overcrowding the .fT-operation is not shown explicitly. 

A formula for P = P 4 (p2), where P 4 represents P 4a or P 4 b, can be readily obtained 
combining Eqs. (12Uj), (1351). and (T4T|): 



M 



4,P 4 (P 2 ) 



dtp( 2 )(t) / (dy) 



P + C Nq + Zq + N l + Z l 



U^Vy 



U 2 V 2 
y y 



U^Vy 



-^\-V y /y a )) 



+ 



N + Z / dU^(-V y /y a ) 
UlVy \ dv y 



(42) 
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where V and Vq are defined by Eq. fl2Tl) . 

In the case P = Pq, where Pq represents one of Pqa, Pqh, we obtain a formula of the 
form 



E + C N + Z N x + Z x 



U 2 V 2 
y y 



-^\-V y /y a )) 



N + Z fdU^(-V y /y a ) 



(43) 



where H^(-V y /y a ) is given by Eq. fl23|) and its derivative can be written in the form 



dU^(-V y /y a ) 



Va 



(dz) 



2-Dp G Bp (\ 
U 2 V 3 U 2 \ V 2 



2Gx 



+ 



C f2x 2Gx< 



- 



Di G B x ( 1 



+ — — + — - + 



Gx\ 



U 3 V 2 U 3 \V V 2 J U A V 



v 3 J u 2 \v 2 



-Vy/ya 



+ 



V 3 



(44) 



As usual the i^-operation and .R-subtraction are assumed implicitly. The i^-operation 
removes UV-divergent parts of divergent subdiagrams from M^p 6a , a = A, B, H , etc. 
Diagrams P§ c and P§d contain fourth-order self-energy subdiagram so that they require 
^-subtraction in addition to i^-operation. The resulting finite quantities are denoted as 
AM 4i p 6o , etc. In order to obtain the standard result renormalized on the mass-shell, further 
subtraction of UV-finite remainder must be carried out by the residual renormalization. 



V. SET 11(C) 

For the Set 11(c) it is convenient to treat the renormalization of UV divergences arising 
from two photons forming M^ a and and the renormalization of the vacuum-polarization 
function separately. The first step can be written as 

A mdMs) [Set n(c)] = J- AM^ } - AB 2 M^ 2) - AB^M 2 , (45) 

i=a,b 

where li,l 2 , and Z3 denote the open lepton line, outer lepton loop, and inner lepton loop, 
respectively. The superscript l\ is suppressed in M 2 and AB 2 since they are independent of 
the lepton mass. AB 2 is the finite part of the second-order renormalization constant defined 
by AB 2 = Bf + Lf. (See Eqs. (|T2|) and (|T4T) .) The vacuum-polarization function Pa{p 2 ) 
is fully renormalized whose divergence structure can be readily found by the i^-operation. 
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This leads to the second step: 

AM iX ( l = E A <SL - ZAB^AM^W, for i = a, 6, (46) 
M £5£, = E M SSl " 2A^ 2 3) M^f , (47) 



A5 £S = E A 4S 2 ) " ZABgf'AB^. (48) 



A. Numerical results: (eee) case 

The contribution of Set 11(c) to the electron g — 2 for the case (hhh) = (eee), where e 
denotes electron, has been evaluated from Eq. fj42l) by three different methods: 



(a) A straightforward extension of the method developed in [291 ], 



hm gencodevpA 



25|, and 



34|. Preliminary eval- 



(b) Method based on the automatic code generating algorit 

(c) Use of an exact spectral function of IT*- 4,2 ) given in Ref. z 
All calculations are carried out by the integration routine VEGAS 
uations of the integral fj42|) by the methods (a) and (b) gave results consistent with each 
other within the numerical uncertainty estimated by VEGAS, proving that both programs 
are bug-free. (Actually, both methods (a) and (b) use only iT-operation since .R-subtraction 
is not needed in this case.) We therefore list only the results of method (b) in the first four 
data lines of Table HI The values of auxiliary functions AB 2 ,p 4{P2) and M 2i p 4(P2) are listed in 
Table ED 

Substituting the numerical results of the integrals listed in Tables HI and HT1 in (f45l) . we 
obtain 

4 10) [Set II(c) (eee) ] = -0.116 489 (32). (49) 
We also calculated the contribution of Set 11(c) using the exact spectral function of II^ 4 ' 2 ) 



28|. In this case AJlfl , AM 4 ( fp } , M 2 ( ee p e) , and AB { 2 ee p e) in the right-hand side of 
Eq. fT4B^ can be directly evaluated using the exact spectral function. The results are listed 
in the last two lines of Table |T] and in the fourth lines of Table [TTJ The value obtained using 
the numbers in Table [I] and Table [Til is 

4 10) [Set H(c) (eee) : spectral function] = -0.116 447 (34) , (50) 
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TABLE I: Contributions of diagrams of Set 11(c) to a e for (hfah) = (eee). The superscript (eee) 
is omitted for simplicity. The multiplicity rip is the number of vertex diagrams represented by the 
integral and is incorporated in the numerical value. The top four lines are obtained by constructing 
the sixth-order vacuum-polarization function in terms of Feynman parameters. The bottom two 
lines are obtained by using the exact spectral function of the sixth order. All integrals are evaluated 
in double precision. 



Integral 


np 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 


AM 4a ,p 4a(P2) 


6 


0.028 927 (21) 


1 x 10 7 , 1 x 10 8 


50, 50 


AM *a,P 4b(P2) 


12 


0.004 521 (11) 


1 x 10 7 , 1 x 10 s 


50, 50 


AM 4b ,p 4a(P2) 


6 


-0.110 617 (16) 


1 x 10 7 , 1 x 10 s 


50, 50 


AM *b,P 4b(P2) 


12 


-0.020 212 (9) 


1 x 10 7 , 1 x 10 8 


50, 50 


Am ^,p 4(P2) 


18 


0.028 425 (28) 


1 x 10 7 , 1 x 10 s 


50, 100 


AM *b,P i(P2) 


18 


-0.112 236 (20) 


1 x 10 7 , 1 x 10 s 


50, 100 



TABLE II: Finite renormalization terms of Set 11(c) for the case (eee). All integrals are evaluated 
in double precision. The quantities in the fourth line are obtained by using the exact spectral 
function of the sixth order. The multiplicity of the integral is incorporated in the numerical value. 



Integral 



Value(Error) 



M 2 ,p 2 
M 2 p„ , a , 

■ 4a(P 2 ) 

M 2 ,p 



Ab(P 2 ) 



Mo 



r 4(P 2 ) 

AM 4fe ,p 2 



0.015 687 421 • 
0.011 403 (1) 
0.001 717 (1) 

0.011 131 (1) 

-0.146 343 (35) 



Integral 



Value (Error) 



2,P 2 



2,-R 



4a(P 2 ) 



AB 2i p 46(P2) 



AS 2,P 4 (P 2 ) 

AM 4a ,p 2 



0.063 399 266 • 
0.047 836 (1) 
0.008 783 (1) 

0.048 577 (5) 

0.039 642 (42) 



which is in good agreement with f|49|) . This shows that GENCODEVP./V works correctly for 
N = 4. 
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TABLE III: Contributions of diagrams of Set 11(c) containing one electron loop and one muon 
loop to a e . The superscript (erne) denotes a diagram in which the outer loop is muon loop. The 
superscript (eem) denotes a diagram in which the inner loop is muon loop. The multiplicity of the 
diagram np is included in the numerical results. All integrals are evaluated in double precision. 



Integral 


np 


Value (Error) 
including np 


Sampling per 
iteration 


No. of 
iterations 


AM ( / m p e) 
AM { d em p e) 

AM< em ;» 

4b,P 4i)( p 2) 


6 
12 
6 

12 


0.000 017 57 (22) 
0.000 013 37 (31) 
-0.000 172 92 (13) 
-0.000 134 62 (18) 


1 x 10 7 
1 x 10 7 
1 x 10 7 
1 x 10 7 


50 
50 
50 
50 


4a >^4a(P 2 ) 

AM^f 
AM^ ] 
AM^ 


6 

12 

6 

12 


0.000 014 66 (11) 
0.000 001 73 (3) 
-0.000 094 42 (8) 
-0.000 001 11 (2) 


1 x 10 7 
1 x 10 7 
1 x 10 7 
1 x 10 7 


50 
50 
50 
50 



B. Numerical results: (erne), (eee), etc. 

Diagrams of Set 11(c) contain two closed lepton loops, one within the other. We obtain 
mass-dependent contributions to the electron g — 2 when one or both loops consist of muon 
or tau-lepton. The largest mass-dependent contributions come from the integral (j45p with 
superscripts {erne) and then with (eem). Results of numerical integration are listed in 
Table IIHI The value obtained using the numbers in Tables II III and [TV] are 

4 10) [Set 11(c) (eme) ] = -0.260 86 (45) x 10" 3 , (51) 

and 

4 10) [Set 11(c) (eem) ] = -0.102 63 (14) x lO" 3 . (52) 
Other mass-dependent terms of Set 11(c) are listed in Table IVl 

C. Muon g — 2: (mee) 

The leading contribution to the muon g—2 comes from the case where both loops consist 
of electrons, namely, the (mee) case, where m stands for the muon. The value obtained 
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TABLE IV: Finite renormalization constants needed for the mass-dependent terms (e, m, e) and 
(e, e, m) of Set 11(c). The finite renormalization constants needed for this term but not listed here 
can be found in Table HH All integrals are evaluated in double precision. The multiplicity of the 
integral is incorporated in the numerical value. 



Integral 



Value (Error) 



Integral 



Value (Error) 



AB. 
AB, 



2,P, 



4a(P 2 ) 



(erne) 
2,Pz 



AB 

AB, 



,^4b(P 2 ) 



(eem) 



2,P, 



4a(P 2 ) 



(eem) 



2,P, 



4b(P 2 ) 



AM, 



(em) 



AB. 
AB, 



ia,P 2 
(em) 

2P2 
(me) 

2,P 2 



0.773 326 (53) x 1(T 4 
0.603 167 (106) x 10~ 4 

0.414 245 (42) x 10 -4 
0.005 000 (9) x 10~ 4 

0.020 90 (21) x 10~ 4 
0.094 050 (3) x 10" 4 
1.885 69 (24) 



M, 



2,P 



4a(P 2 ) 



M, 



(erne) 



2,P 



46(P 2 ) 



,1/ 



(eem) 



2,P 



4a(P 2 ) 



M, 



(eem) 



2,P 



4b(P 2 ) 



(em) 



M (em) 



0.044 986 (3) x 10" 4 
0.033 465 (6) x 10~ 4 

0.050 072 (5) x 10" 4 
0.000 552 (1) x 10~ 4 

-0.209 84 (12) x 10~ 4 
0.005 197 62 (21) x 10" 4 



TABLE V: Mass-dependent contributions of diagrams of Set 11(c) to the electron g— 2. All integrals 
are evaluated in double precision. 



A 



(10)(eZ 2 Z 3 ) 



A 



(10)(el 2 l 3 ) 



(e, m, m) 
(e,t,e) 
(e,e,t) 
(e,t,t) 



-0.167 65 (28) x 10~ 4 
-0.287 97 (58) x 10~ 5 
-0.988 9 (20) x 10" 6 
-0.988 4 (21) x 10~ 7 



(e,m,t) 
(e,t,m) 



-0.410 01 (81) x 10~ 6 
-0.784 4 (14) x 10~ 6 



using the numbers in Tables I VI I and IVHI is 

4 10) [Set 11(c) (mee) ] = -3.888 27 (90). 
We checked this result using the exact spectral function: 

4 10) [Set II(c) (mee) : spectral function] = -3.887 65 (92). 



(53) 



(54) 



D. Muon g — 2: (mme), (mem) and others 

The next-to-leading order contribution arises when the inner and outer loops consist of 
electron and muon, respectively. We found 

4 10) [Set 11(c) (mme) ] = -1.345 98 (36). (55) 
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TABLE VI: Contributions to the muon g — 2 from Set 11(c) diagrams involving closed electron 
and/or muon loops. The multiplicity of the diagram np is included in the numerical results. All 
integrals are evaluated in double precision. 



Integral 


Up 


Value (Error) 
including np 


Sampling per 
iteration 


No. of 
iterations 


AMjr e p e) 

AM (mee) 
^ IV ha,P 4b(P2) 

AM (mee) 
AM (mee) 


6 
12 
6 

12 


0.684 47 (37) 
2.071 36 (55) 
0.025 50 (33) 
-4.320 77 (51) 


1 x 10 8 , 1 x 10 9 
1 x 10 8 , 1 x 10 9 
1 x 10 8 , 1 x 10 9 
1 x 10 8 , 1 x 10 9 


50, 40 
50, 50 
50, 40 
50, 50 



AM (mme) 

AM (mme) 
^ M 4a,P 4b(P2) 

AM 4b,P 4o(P2) 

AM 4b,P 4b(P2) 

AM (mem) 

A M ( mem ) 
^ IV ha,P 4b{P2) 

A ]\A mem ) 
^ M 4b,P 4a{P2) 

^ M 4b,P 4b(P2) 



6 

12 

6 

12 

6 

12 
6 

12 



0.302 14 (11) 
0.251 36 (19) 
-0.957 39 (9) 
-0.931 56 (15) 

0.049 57 (10) 
0.001 64 (5) 
-0.141 38 (8) 
-0.010 38 (4) 



1 x 10 8 
1 x 10 8 
1 x 10 8 
1 x 10 8 



1 x 10 9 
1 x 10 9 
1 x 10 9 
1 x 10 9 



1 x 10 8 
1 x 10 8 
1 x 10 s 
1 x 10 s 



50, 10 
50, 10 
50, 10 
50, 10 

50 
50 
50 
50 



When the inner and outer loops consists of muon and electron, respectively, the contribution 
is found to be smaller: 

4 10) [Set II(c) (mem) ] = -0.151 50 (15). (56) 
The contributions involving tau-lepton loops are summarized in Table IVIIII 

VI. SET 11(D) 

Following the same consideration leading to Eq. ( I4"5j) of Sec. |V]we obtain 

b 

4 10) [Set II(d) (Zl/2) ] = AM ^)e - AB 2 M 2 h pf - A5^ 6 2) M 2 , (57) 

i=a 

where I2 designates the loop lepton. The vacuum-polarization function P$ is fully renormal- 
ized whose divergence structure can be readily found by the i^-operation. The renormal- 



18 



TABLE VII: Finite renormalization terms of Set 11(c) for the muon anomaly a„. All integrals 
are evaluated in double precision. The multiplicity of the integral is incorporated in the numerical 
value. 



Integral Value(Error) Integral Value(Error) 



AB^ m p ee) 
ABi m p ee) 

> 4b(P 2 ) 
2,P 4 (P 2 ) 


0.655 71 (11) 
2.279 41 (17) 
2.695 12 (64) 


jy-imee) 

2 , P 4a(P 2 ) 

M 2 ( T e) 

> 4b(P 2 ) 

jy-imee) 
2,P 4 (P 2 ) 


0.597 44 1 (48) 
0.982 06 6 (70) 
1.440 46 (28) 


AB (mme) 

2 , P 4a(P 2 ) 

a gimme) 

2 , P 4b(P 2 ) 


0.417 691 (10) 
0.404 336 (22) 


M (mme) 

2 , P 4a(P 2 ) 

2 , P 4b{P 2 ) 


0.121 908 (3) 
0.099 237 (5) 


2 , P 4a(P 2 ) 

^£,(mem) 

2 , P 4b(P 2 ) 


0.065 066 (6) 
0.004 533 (3) 


2 , P 4a(P 2 ) 
2 > P 4b(P 2 ) 


0.021 016 (2) 
0.000 586 (1) 


amW 

lSM 4a,P 2 
ad (me) 
^ n 2,P 2 


1.725 62 (49) 
1.885 732 (16) 


AM (me) 

M (me) 


-2.354 33 (45) 
1.094 258 282 7 (98) 



TABLE VIII: Contributions to the muon g— 2 from Set 11(c) diagrams involving tau-lepton loops. 
All integrals are evaluated in double precision. 





A (10)(ml 2 l 3 ) 




A (10)(ml 2 l 3 ) 


(m, m, t) 


-0.004 325 (49) 


(m,e,t) 


-0.004 734 1 (55) 


(m, t, m) 


-0.010 519 (13) 


(m,t,e) 


-0.036 066 (51) 


(rn,t,t) 


-0.001 504 1 (19) 







ization formula for P 6 takes different forms depending on whether one follows the original 
.fT-operation prescription 29j or the i^-operation plus .R-subtraction method 19] . In the 
first approach the UV-finite part of fourth-order mass-renormalization term is not sub- 
tracted when Pqc and Pqd are constructed. In the second approach we subtract the mass- 
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renormalization term completely, including the finite part ASm^, which leads to 

A <£ = E A <£ 

- AAB 2 AM^X ] ~ S(AB 2 ) 2 AM^ - 2ALB 4 AM^, for i = a, b , 

<k ] = E M S 

I3=A 

- AAB 2 M ( 2 h pf - 3(A5 2 ) 2 M 2 ( ^ 2 2) - 2ALB i M { 2 h p l f , 
*B$g> = E ABgg 

- AAB 2 AB { 2 h pf - 3(AB 2 ) 2 AB^pf - 2ALB 4 AB { 2 ff . (58) 

The quantities in the right-hand-side of (J58"]l . AM^, A£? 2 ,p 6/3 , and M 2i p 6)9 are defined by 
the A"-operation and i?-subtraction. ALB 4 is the sum of the finite parts of the fourth-order 
vertex-renormalization constant AL4 and wave-function renormalization constant AB 4 . See 
Refs. QBE I for the exact definition. Note that terms like AM^p 60 include the multiplicity 
/If of Feynman diagrams that contribute to them. 

A. Numerical results: (ee) case 

Preliminary calculations of the (ee) case based on Methods (a) and (b) described in 
Sec. IV Al are consistent with each other within the uncertainty estimated by VEGAS. There- 
fore we list only the results of method (b) in Table IIXI From this table and Table |X] we 
obtain 

4 10) [Set 11(d) (ee) ] = -0.243 00 (29) . (59) 



The Pade-approximated vacuum-po 
fermion loop has been obtained in Ref. 



arization function of the sixth-order with a single 
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31] . This method gives both imaginary and real 



parts of the vacuum-polarization function. We use here only its imaginary part to calculate 
its effect on the anomaly. Numerical results of integration are summarized in Table. IXHI 
Substituting them into Eq. (137)1 . we obtained 



4 10) [Set 11(d) (ee> : Pade] = -0.243 06 (45) 



\( ee ) . 



(60) 
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TABLE IX: Contributions of diagrams of Set 11(d), (ee) case, np is the number of Feynman 
diagrams represented by the integral. The fourth-order mass-renormalization is completed by R- 
subtraction within the numerical programs of AM 4a) p 6C , AM4 ai p w , AM^p 6C , and AM^p 6D . All 
integrals are evaluated in double precision. 



Integral 


np 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 





12 


0.112 


990 


(116) 


1 X 


10 8 , 


1 


x 


10 9 


50, 


100 


AMdn P 


6 


0.072 


919 


(72) 


1 X 


10 8 , 


1 


x 


10 9 


50, 


100 


AMdn P 


12 


0.044 


224 


(87) 


1 X 


10 8 , 


1 


x 


10 9 


50, 


100 


AMa„ p a „ 

'— »-"'- t 4a,-f'6.D 


12 


-0.088 


822 


(78) 


1 X 


10 8 , 


1 


x 


10 9 


50, 


100 




21 


0.444 


033 


(113) 


1 X 


10 8 , 


1 


X 


10 9 


50, 


100 


AM 4a ,p 6F 


12 


-0.156 


407 


(67) 


1 X 


10 8 , 


1 


X 


10 9 


50, 


100 


AM 4a ,p 6G 


6 


0.094 


162 


(54) 


1 X 


10 8 , 


1 


X 


10 9 


50, 


100 


AM *a,P 6H 


6 


0.060 


989 


(35) 


1 X 


10 8 , 


1 


X 


10 9 


50, 


100 




12 


-0.398 


926 


(66) 


1 X 


10 8 , 


1 


X 


10 9 


50, 


100 


AM *b,p 6B 


6 


-0.253 


369 


(42) 


1 X 


10 8 , 


1 


X 


10 9 


50, 


100 




12 


-0.141 


941 


(51) 


1 X 


10 8 , 


1 


X 


10 9 


50, 


100 


AM 4b ,p 6D 


12 


0.292 


773 


(44) 


1 X 


10 8 , 


1 


X 


10 9 


50, 


100 


AM 4b ,p 6E 


21 


-1.395 


971 


(66) 


1 X 


10 8 , 


1 


X 


10 9 


50, 


100 


AM 4b ,p 6F 


12 


0.570 


363 


(40) 


1 X 


10 8 , 


1 


X 


10 9 


50, 


100 


AM 4fe,P 6G 


6 


-0.232 


467 


(32) 


1 X 


10 8 , 


1 


X 


10 9 


50, 


100 




6 


-0.223 


983 


(22) 


1 X 


10 8 , 


1 


X 


10 9 


50, 


100 



which is in good agreement with ( 1591) . This provides another support for the validity of 

GENCODEVPA. 

B. Numerical results: mass-dependent terms (em) and (et) 

The value of the mass-dependent term (em) obtained using the numbers listed in Table IXIl 

is 

4 10) [Set 11(d) (em) ] = -0.981 7 (42) x 10" 4 . (61) 

As a check we evaluated the same quantity using the Pade-approximated vacuum- 
polarization function of sixth-order. The results are listed in Table IXHI From these values 
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TABLE X: Finite renormalization terms necessary for the cases (ee) and (em) of Set 11(d). For 
simplicity the superscript (ee) is omitted. The values M 2 ^p ec and M 2 ^p eD are different from those 



in Table. I of Ref. 



351 ] . The former is constructed with the JC-operation and .R-subtraction, while 



the latter is with the ET-operation only. All integrals are evaluated in double precision. 



Integral 



Value(Error) 



Integral 

AM 46i p 4 
AM 4b:P2 
AB 2>Pf . A 
AB 2 ,p 6B 

AB 2,P 6C 
AB 2,P 6E 

A-B 2 ,P 6H 

A5 2 ,p 4 
AB 2t p 2 
AB 2 



Value(Error) 



AM 4o ,p 4 
AM 4a ,p 2 

M2,P 6A 
M2,P 6B 
M2,P 6C 
M2,P 6D 
M2,P 6E 
M2,P 6F 

M2,P 6G 
M2,P 6H 

M 2 ,p 4 

M 2 ,p 2 
ALB A 
M 2 



AM (em) 
AM (em) 

<p£ 

mH 

IV1 2,P 6E 
M (em) 

2,P 6 G 
2,P6iT 

mW 

JW 2,P 4 

M (em) 

JW 2,P 2 



0.131 298 (8) 

0.039 642 (42) 

0.044 446 7 (22) 

0.028 593 9 (14) 

0.017 717 3 (19) 
-0.035 167 (16) 

0.179 333 2 (21) 
-0.062 003 2 (12) 

0.038 879 (10) 

0.023 674 9 (8) 

0.052 870 652- •• 

0.015 687 421 •• • 

0.027 930 (27) 

0.5 

0.075 96 (78) x 10~ 4 
0.020 90 (21) x 10" 4 
0.159 949 (32) x 10~ 5 
0.102 323 (21) x 10~ 5 
0.071 033 (29) x 10~ 5 

-0.131 368 (26) x 10~ 5 
0.682 404 (32) x 10~ 5 

-0.211 985 (19) x 10~ 5 
0.174 650 (16) x 10~ 5 
0.082 748 (12) x 10~ 5 
0.197 298 (5) x 10" 5 
0.051 974 (1) x 10~ 5 



AM 



(em) 



4b, P 4 



AM 



(em) 



AB. 
AB, 
AB, 



4b,P 2 
(em) 

2,P 6 A 
(em) 

2,P 6 s 
(em) 

'2,P 6C 

ABi e p' 
AB^f 

AB^l 
ABif 



-0.420 295 (8) 

-0.146 343 (35) 
0.173 609 1 (96) 
0.110 466 1 (56) 
0.062 134 7 (74) 

-0.127 657 6 (63) 
0.610 385 8 (84) 

-0.247 658 3 (54) 
0.104 070 2 (44) 
0.097 056 7 (29) 
0.183 666 8 (18) 
0.063 399 266 ■■ • 
0.75 



-0.757 35 (41) x 10~ 4 

-0.209 84 (12) x 10" 4 
0.281 371 (61) x 10~ 4 
0.180 003 (40) x 10~ 4 
0.119 848 (52) x 10~ 4 

-0.226 651 (46) x 10~ 4 
1.162 709 (59) x 10~ 4 

-0.379 345 (35) x 10~ 4 
0.279 935 (31) x 10~ 4 
0.147 672 (22) x 10~ 4 
0.338 738 (12) x 10~ 4 
0.094 050 (3) x 10~ 4 



we obtain 

4 10) [Set 11(d) (em) : Pade] = -0.991 53 (61) x 1(T 4 . 
We also evaluated the mass- dependent term (et) in Pade approximation: 

4 10) [Set 11(d) (et) : Pade] = -0.542 7 (18) x 10~ 6 . 



(62) 



(63) 
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TABLE XI: Contributions to a e from diagrams of Set 11(d) containing a muon loop. The super- 
script (em) signifies that the diagrams contain muon loop in the electron g — 2. np is the number 
of Feynman diagrams represented by the integral. All integrals are evaluated in double precision. 
The fourth-order mass renormalization is completed within the numerical programs of AM^p 6C , 
AM 4a ,p 6D , AM 4fejP6C , and AM 4bi p 6D . 



Integral 


np 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 



(em) 



AM, 



4a,P 6A 



AM, 



(em) 



AM, 



(em) 



4a,P 6C 
4a,P 60 

A <1 

a <pL 
A <1 

A <1 

A <t 

Jem) 

A <1 

r (em) 



am! 



AM 



4b,P fi 



12 


0.000 


006 


24 


(13) 


1 X 


10 7 , 


1 


X 


10 8 


50, 


50 


6 


0.000 


004 


94 


(9) 


1 X 


10 7 , 


1 


X 


10 8 


50, 


50 


12 


0.000 


003 


79 


(11) 


1 X 


10 7 , 


1 


X 


10 s 


50, 


50 


12 


0.000 


005 


23 


(10) 


1 X 


10 7 , 


1 


X 


10 s 


50, 


50 


24 


0.000 


026 


34 


(15) 


1 X 


10 7 , 


1 


X 


10 s 


50, 


50 


12 


0.000 


008 


38 


(8) 


1 X 


10 7 , 


1 


X 


10 8 


50, 


50 


6 


0.000 


007 


54 


(6) 


1 X 


io 7 , 


1 


X 


10 8 


50, 


50 


6 


0.000 


003 


34 


(5) 


1 X 


io 7 , 


1 


X 


10 8 


50, 


50 


12 


0.000 


063 


78 


(4) 


1 X 


io 7 , 


1 


X 


10 8 


50, 


50 


6 


0.000 


040 


16 


(2) 


1 X 


10 7 , 


1 


X 


10 8 


50, 


50 


12 


0.000 


027 


79 


(3) 


1 X 


10 7 , 


1 


X 


10 8 


50, 


50 


12 


0.000 


051 


62 


(3) 


1 X 


10 7 , 


1 


X 


10 s 


50, 


50 


24 


0.000 


260 


81 


(4) 


1 X 


10 7 , 


1 


X 


10 8 


50, 


50 


12 


0.000 


085 


61 


(2) 


1 X 


io 7 , 


1 


X 


10 8 


50, 


50 


6 


0.000 


063 


72 


(2) 


1 X 


10 7 , 


1 


X 


10 8 


50, 


50 


6 


0.000 


033 


95 


(1) 


1 X 


io 7 , 


1 


X 


10 8 


50, 


50 



We have not evaluated this term directly. But it will be of the same order as Eq. (|63|) and 
thus negligible numerically. 

C. Muon g — 2: (me) 

The leading contribution to the muon g—2 comes from the case (me). The value obtained 
using the numbers in Tables IXHIl IXl and IXVI is 

4 10) [Set 11(d) (me) ] = 0.497 2 (65). (64) 

This is in fair agreement with the value obtained using the Pade approximant 

4 10) [Set 11(d) (me) : Pade] = 0.504 8 (75). (65) 
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TABLE XII: Contributions of diagrams of Set 11(d) whose VP function is Pq. Quantities on the 
right-hand-side of Eq. (ESty are calculated from Tables US H EH IXTO1 IXTVl and |XVJ The same 
quantities are also calculated by the Pade approximant method and listed with the subscript Pade 
below the corresponding integrals. M^pQ are actually the anomaly contributions of the eighth- 
order diagrams of Group 1(d). Their values for the (ee) and (me) cases are consistent with Eq. (29) 

n 

of Ref. [35f] and Eq. (34) of Ref. [a], respectively. The (em), (et), and (mt) cases are newly evaluated 



in this paper. All integrals using the Pade approximant are evaluated in quadruple precision. 



Integral 



Value(Error) 



Integral 



Value (Error) 



AM 



4a,P 6 



AM 



(ee) 



AB. 
AB, 



4a, Pq Pade 

(ee) 



2,-Pe 
(ee) 



2,P 6 Pade 



0.119 08 (25) 
0.118 95 (35) 
0.120 879 (20) 
0.120 862 (39) 



AM 



4fc,P 6 
^ iW 4&,P 6 Pade 

M (ee) 

M (ee) 
1V1 2,P 6 Pade 



-0.264 51 (41) 
-0.264 43 (27) 
0.049 514 (5) 
0.049 520 (4) 



AM, 



(em) 



4a,P 6 



AM, 



(em) 



AS 2,P 6 

abH 



4a, P 6 Pade 
(em) 



2,P 6 Pade 



0.091 5 (39) x 10~ 4 
0.091 96 (53) x 10~ 4 
0.385 36 (13) x 10" 4 
0.385 367 (72) x 10~ 4 



^ iW 46,P 6 Pade 

(em) 

2,P 6 
(em) 

2,P 6 Pade 



M 
M 



-0.862 (16) x 10~ 4 
-0.872 25 (31) x 10~ 4 
0.024 725 (7) x 10~ 4 
0.024 727 (4) x 10~ 4 



AM, 



(et) 



AB. 



4a. P 6 Pade 

(et) 

2,P 6 Pade 



0.039 89 (154) x 10" 
0.210 278 (39) x 10" 



LS1V1 4b,P (i Pade 

Aet) 

J 2,P 6 Pade 



-0.470 90 (74) x 10" 
0.008 744 (1) x 10" 



AM, 



(me) 



4a,P 6 

^■ Ivl 4a,P 6 Pade 

AB (me) 

^ A - D 2,P 6 

£A - D 2,P 6 Pade 



-0.469 4 (50) 

-0.459 1 (55) 
-0.394 72 (82) 
-0.395 25 (78) 



AM, 



(me) 



46,P 6 
^ AiW 4fe,P 6 Pade 
1V1 2,P 6 
1V1 2,P 6 Pade 



0.596 9 (42) 
0.593 5 (51) 
-0.229 82 (36) 
-0.230 23 (32) 



AM 



(mt) 



4a,P 6 



AM, 



(mt) 



AB. 
AB. 



4a, P 6 Pade 
(mt) 



2,P 6 
(mt) 

2,P 6 Pade 



0.001 066 (15) 
0.001 065 (12) 
0.003 020 9 (10) 
0.003 021 4 (6) 



AM, 



(mt) 



4fe,P 6 

^ iW 46,P 6 Pade 
M (mt) 
iW 2,P 6 
M (mt) 
iW 2,P 6 Pade 



-0.006 953 (10) 
-0.006 953 (8) 
0.000 367 7 (2) 
0.000 367 7 (1) 



A crude evaluation of the contribution of tau-lepton loop gives 



4 10) [Set 11(d) 



(mf)-i 



-0.007 673 (18), 



(66) 
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TABLE XIII: Leading contributions of diagrams of Set 11(d) of Fig. Q] to the muon g— 2. np is the 
number of Feynman diagrams represented by the integral. All integrals are evaluated in double 
precision. The fourth-order mass-renormalization is completed within the numerical programs of 
AM4 aj p 6C ,, AM4a t p aD , AM^Pgp, and AM^)P 6D . Last two columns list initial sampling and its 
iteration, followed by increased sampling and its iteration. 



Integral 


np 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 



AM (me) 
A <t 

AM (me) 
AM (me) 

A <1 

am'™' 

r(me) 

Ub,p 6F 
AM (me) 

^ 1V1 4b,P 6G 

AM (me) 

^ IVI 4b,P 6H 



AM; 



12 


10.693 


5 (20) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


185 


6 


5.642 


1(12) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


170 


12 


4.305 


(12) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


170 


12 


-6.097 


7(12) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


170 


24 


-4.365 


9 (18) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


170 


12 


-6.145 


99 (94) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


170 


6 


-0.397 


7(14) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


170 


6 


5.048 


76 (52) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


100 


12 


-18.959 


1(16) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


170 


6 


-9.815 


33 (96) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


150 


12 


-7.998 


86 (93) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


150 


12 


11.248 


41 (81) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


150 


24 


15.905 


(12) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


170 


12 


7.177 


72 (74) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


150 


6 


-0.501 


22 (94) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


160 


6 


-8.025 


32 (40) 


1 


X 


10 8 , 


1 


X 


10 9 , 


1 


X 


10 10 


50, 


100, 


100 



while the same contribution obtained using the Pade approximant is 

4 10) [Set 11(d) {mt) : Pade] = -0.007 674 (15), 



(67) 



VII. SUMMARY AND DISCUSSION 

The total contribution to a e from the set 11(c) is the sum of Eqs. (l4"9|) . ( l5"Tj) . ( 152|) . and 
other terms listed in Table |Vj 

a( 10 )[Set II(c):aU] = -0.116 874 (43) (^) 5 . (68) 

Contributions of tau-lepton loop listed in Table rVTTTI are less than the uncertainty of Eq. fl68|) . 
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TABLE XIV: Contribution from diagrams of Set 11(d) (to, t) of Fig. [T] to the muon g— 2. np is the 
number of Feynman diagrams represented by the integral. All integrals are evaluated in double 
precision. Last two columns list initial sampling and its iteration, followed by increased sampling 
and its iteration. 



Integral 


np 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 



AM 



{rat) 



4a,P 6A 



AM, 



(rat) 



4a,P 6s 



AM, 



(rat) 



AM) 



Jrat) 
L Aa,P 6F 

AM (mt) 

^ IV1 4b,P 6C 

AM (mt) 

^ 1V1 4b,P 6D 



12 


0.001 


746 1 (22) 


1 


x 


10 s , 


1 


x 


10 9 


50, 


100 


6 


0.000 


476 4 (14) 


1 


x 


10 8 , 


1 


x 


10 9 


50, 


100 


12 


0.000 


323 9 (18) 


1 


x 


10 8 , 


1 


x 


10 9 


50, 


100 


12 


-0.001 


605 3 (17) 


1 


x 


10 8 , 


1 


x 


10 9 


50, 


100 


24 


003 


112 4 


1 




10 8 


1 




10 9 


50 


1 00 


12 


-0.002 


000 (13) 


1 


X 


10 8 , 


1 


X 


10 9 


50, 


100 


6 


0.001 


766 5 (12) 


1 


X 


10 8 , 


1 


X 


10 9 


50, 


100 


6 


0.000 


390 55 (72) 


1 


X 


10 8 , 


1 


X 


10 9 


50, 


100 


12 


-0.006 


549 50 (87) 


1 


X 


10 8 , 


1 


X 


10 9 


50, 


100 


6 


-0.004 


550 21 (57) 


1 


X 


10 8 , 


1 


X 


10 9 


50, 


100 


12 


-0.002 


299 84 (74) 


1 


X 


10 8 , 


1 


X 


10 9 


50, 


100 


12 


0.004 


413 20 (67) 


1 


X 


10 8 , 


1 


X 


10 9 


50, 


100 


24 


-0.022 


437 50 (93) 


1 


X 


10 8 , 


1 


X 


10 9 


50, 


100 


12 


0.008 


566 26 (53) 


1 


X 


10 8 , 


1 


X 


10 9 


50, 


100 


6 


-0.005 


164 07 (46) 


1 


X 


10 8 , 


1 


X 


10 9 


50, 


100 


6 


-0.003 


940 06 (31) 


1 


X 


10 8 , 


1 


X 


10 9 


50, 


100 



The total contribution to a e from the Set 11(d) is the sum of Eq. (|59j) and Eq. (IBTjl : 

a( 10) [Set II(d):all] = -0.243 10 (29) (^) 5 . (69) 

The contribution of the tau-lepton loop is within the error bars of Eq. (|69|) and is completely 
negligible at present. 

The total contribution of Set 11(c) to the muon g — 2 involving electron, muon, and tau- 
lepton loops is the sum of Eq. ( )49|) . Eq. ( 1531) . Eq. ( 1551) . and Eq. ( 1561) . and values listed in 
Table. M 

a^ 10) [Set II(c):all] = -5.559 4 (11) (-) 5 . (70) 
The total contribution of Set 11(d) to the muon g — 2 involving electron, muon, and 
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Value (Error) 


— Z.40D 


fin 
DU 


(,12J 


— Z.O04 


33 


(46) 


11.500 


82 


(47) 


6.101 


83 


(27) 


4.616 


02 


(24) 


-6.701 


64 


(22) 


-4.076 


72 


(36) 


-6.535 


49 


(20) 


-0.039 


85 


(28) 


5.345 


15 


(8) 


2.439 


109 (53) 


1.885 


733 (16) 



TABLE XV: Finite renormalization terms needed for Set 11(d) of the muon g— 2. Bothfme) and 

(mt) cases are listed. M^p^ with f3 = A,B, E, F, G, H are consistent with those in Ref. [5J. Mop 6i3 

with /3 = C, D in this table incorporated the i?-subtraction and differs from those in Ref. 5]. 

Integral Value(Error) Integral 

AM^ P { 2.047 84 (23) AM]™j 

AM^ 1.725 62 (49) AM i™$ 2 

5.676 49 (22) AB 2\Pe A 

M^J 3.058 13 (13) AB 2 m ptB 

M^pJ c 2.194 66 (12) as £f? c 

M 2 m pf D -3.224 25 (10) A - B S ) D 

M 2 ( ^ -0.073 76 (17) AB< %\pJ 

M 2^ F -4.064 09 (9) AB tp% 

-0.246 97 (12) as £Sg 
M 2 ( ^ 2.838 67 (4) 

M 2 ( ^ 1.493 671 581 (8) ^tp} 

M^pf 1.094 258 282 7 (98) AB £r? 

AM^ 4 0.000 903 9 (47) AM i™p 4 -0.006 571 6 (31) 

AM^ 2 0.000 247 8 (15) AM 4™% -0.001 888 8 (10) 

0.000 239 71 (1) AB< Cp] a °- 002 436 62 ( 13 ) 

M 2 ( ™* } s 0.000 153 39 (1) AR SSL °- 001 558 94 ( 8 ) 

0.000 106 19 (1) AB{ 2 m p ] c °- 001 006 23 ( n ) 

-0.000 196 75 (1) AB 2\p d -0.001 934 79 (9) 

M 9 ( ™* } 0.001 021 35 (1) AB^ 0.009 826 56 (12) 

M 2 ( ™* } f -0.000 318 06 (1) AB 2\Pg F -0.003 326 19 (7) 

M 2 ( ™* } g 0.000 260 44 (1) AR Ste °- 002 250 09 ( 6 ) 

0.000 124 05 (1) AB 2!ps H °- 001 292 36 ( 4 ) 

M^p] 0.000 295 508 (21) AB t}\ °- 002 880 01 ( 31 ) 

M { ™p] 0.000 078 067 4 (31) AB { ™p] °- 000 831 107 ( 75 ) 

tau-lepton loops is the sum of Eq. ( 1591) . Eq. ( 1(541) . and Eq. ( 1(561) : 

aj^fSet II(d):all] = 0.246 5 (65) (-) 5 . (71) 



The sum of the electron-loop contribution to the muon g— 2 from the diagrams Set 11(c) 
and Set 11(d) is the sum of Eqs.([53]), (j55]b (|55|), and (jM]). We find 

4 10) (m M /m e )[Set II(c+d) (me) ] = -4.888 6 (65), (72) 
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which is less than 1 % of the leading contribution from the diagrams of Set VI (a) that contain 



light-by-light-scattering subdiagrams and vacuum-polarization subdiagrams 
the new contribution does not alter the previous estimates: 



A 



(10), 



t /m e ) [estimate : Ref. 



663 (20) 



A 



(10), 



m M /m e ) [estimate : Ref.^g]] = 643 (20) . 



Hence, 

(73) 
(74) 
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